A direct measurement of the Bose-Einstein Condensation universality class in 
NiCl 2 -4SC(NH 2 )2 at ultra-low temperatures 
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In this work, we demonstrate field-induced Bose-Einstein condensation (BEC) in the organic 
compound NiCl2-4SC(NH2)2 using AC susceptibility measurements down to 1 mK. The Ni S=l 
spins exhibit 3D XY antiferromagnetism between a lower critical field H c i ~ 2 T and a upper 
critical field H C 2 ~ 12 T. The results show a power-law temperature dependence of the phase 
transition line H cl (T) - H cl (0) = aT a with a = 1.47±0.10 and H cl (0) = 2.053 T, consistent with 
the 3D Bose-Einstein Condensation universality class. An abnormal change was found in the phase 
boundary near H C 2 at approximately 150 mK. 



The idea of Bose-Einstein condensation (BEC) occur- 
ring in the spin systems of certain quantum magnets with 
axial symmetry has been explored extensively in the past 
few years This idea was first suggested by Affleck in 
1991 0] and first investigated experimentally in the com- 
pound TlCuCl 3 [3]. Although BEC in TlCuCl 3 was later 
called into question in an oft-overlooked electron spin res- 
onance (ESR) study [4[ , the work on TICUCI3 generated 
a flurry of experimental and theoretical activity in the 
field and many new candidate BE C sy stems have since 
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been proposed [a, la IH, |8|, |9j, 

In these compounds, the axial symmetry of the spins 
allows XY antiferromagnetic order to occur over certain 
ranges of temperature and magnetic field. Near the criti- 
cal magnetic fields where the long-range order is induced 
or suppressed, the spin system can be mapped onto a 
system of dilute hard-core bosons on a lattice and the 
field-induced quantum phase transition at the bound- 
ary of the long-range ordered state can be modeled as 
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a Bose-Einstein condensation 
caveat is that in quantum magnets the boson number 
is proportional to the longitudinal magnetization, and 
thus is conserved in equilibrium rather than strictly on 
all time scales. Therefore, in quantum magnets only 
the thermodynamic properties of the system should fol- 
low the predictions of BEC theory and nonequilibrium 
effects such as supercurrents will not occur. Neverthe- 
less, these compounds provide an important test of BEC 
phase transition in the thermodynamic limit. 

A key experimental signature of BEC is a power- 
law temperature-dependence of the number of condensed 
bosons with an exponent of 3/2. This power-law is the 
low-temperature limit of the boson distribution function. 
In the spin systems this translates to a power-law tem- 
perature dependence of the critical field line H c — H c (0) ex 
T a where a — 3/2 (I9L lioj]. It is important to note 
that this power-law is valid in the limit of very low tem- 
peratures. So far very few quantum magnet BEC candi- 



date have been studied at temperatures well below the 
energy scales for boson interactions, which are given by 
the antiferromagnetic couplings between spins. Further- 
more, if the temperature range at which the power-law fit 
is performed is far from zero temperature, then it is very 
difficult to accurately identify the intercept, and the re- 
sulting power-law exponent that is derived from the fit is 
highly dependent on the value of the intercept used 21 1 . 
One way to circumvent this problem is the windowing 
method in which the intercept and the exponent are de- 
termined more-or-less independently by performing fits 
over different temperature ranges and extrapolating the 
values of the intercept and the exponent to zero temper- 
ature 22]. However, this extrapolation technique from 
higher temperatures does not always correspond to the 
actual low-temperature behavior as was seen in the com- 
pound BaCuSi20g, where at temperatures above 1 K, the 
windowing method yields one exponent, but due to a re- 
duction in dimensionality of the spin system, a different 
exponent can be observed at lower temperatures [Hj]. 

Here we report a direct observation of the 3/2 power 
law exponent of the 3D BEC universality class at ul- 
tra low temperatures in the compound DTN (MCI2- 
4SC(NH2)2)- Our measurements were performed at tem- 
peratures down to 1 mK, which is two orders of mag- 
nitude below the lowest temperature scale for magnetic 
coupling in this system J a — 180 mK, and the lowest 
temperature ever used to investigate BEC in a quantum 
magnet. Thus we do not have to use any extrapolation 
to determine the power-law exponent. 

The organic magnet DTN contains 5=1 Ni 2+ atoms 
that form two interpenetrating tetragonal lattices. At 
zero field, a uniaxial anisotropy D~9K splits the Ni S 
= 1 triplet into a S z — ground state and a S z = ±1 ex- 
cited doublet. The S z = 1 state can be suppressed with 
applied magnetic fields along the tetragonal c-axis via the 
Zeeman effect, thus producing a magnetic ground state 
above H c i = 2.1 T. Antiferromagnetic coupling between 
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the Ni atoms with exchange strength J c = 1.8 K along 
the c-axis and J a = 0.18 K perpendicular to the c-axis 
produces long-range antiferromagnetic order [24]|. The 
long-range order occurs in a dome-shaped region of the 
T — H phase diagram between H c \ , where the magnetic 
ground state is induced, and H c2 where the spins align 
with the applied magnetic field and below T = 1 K 7]. 
The XY symmetry, the magnetic exchange couplings and 
the uniaxial anisotropy D have been identified via in- 
elastic neutron scattering and electron spin resonance 
(ESR). The power-law dependence of the critical field 
line H c \ was previously investigated using specific heat, 
magnetocaloric effect and magnetization measurements 
down to a minimum temperature of 100 mK [7H2H] using 
the windowing method, and a power-law consistent with 
a = 3/2 was determined. The validity of this power-law 
was called into question however due to the fact that it 
is an extrapolation rather a direct measurement in par- 
ticular in light of the recent new results for BaCuSi^Oe 
showing a different exponent at low temperatures. [23| 

We have now measured the power-law temperature de- 
pendence of the critical field H& down to 1 mK using AC 
susceptibility measurements. The experiments were car- 
ried out using a PrNis nuclear refrigerator and a 15 T 
magnet at the High B/T facility of National High Mag- 
netic Field Laboratory. The sample was immersed in 
liquid 3 He in a polycarbonate cell, and thermal contact 
to the refrigerator was assured via sintered silver that 
was an integral part of an assembly of annealed silver 
rods extending form the nuclear refrigerator. The tem- 
perature was calibrated by a 3 He melting-pressure curve 
thermometer mounted in the zero-field region of the mag- 
net, at the top of the nuclear stage. Both AC and DC 
magnetic fields were applied in the direction parallel to 
the c-axis of the single crystal of DTN. An AC signal 
with an amplitude of 0.5 ~ 1 G and a low frequency 
of 275 Hz was generated in a primary coil wound from 
superconducting NbTi wire, thereby avoiding heating at 
ultra-low temperatures. The AC susceptibility Xac was 
measured by sweeping the external DC field at rates be- 
tween 0.1 and^ 0.001 T/min while the temperature was 
fixed at various values between 270 and 1.0 mK. 

Fig. Q] shows three traces of AC susceptibility \ac for 
temperatures of 20, 150, and 600 mK with a relatively 
fast sweep rate of 0.054 T/min. The typical field-induced 
long-range AFM (BEC) transitions appears as steps in 
these traces with H& ~ 2 T and H C 2 ~ 12 T. As the 
temperature is lowered, the steps become sharper and a 
peak develops near H C 2. 

The values of H c \ and H C 2 were determined from 
the peak in the first derivative of the AC susceptibility 
dXac/d-ff , as shown in the inset of Fig. [2j In the critical 
field - temperature phase diagram shown in Fig. [2]^a) and 
Fig. [2fb), the temperature dependence of H c \ and H c2 
have been plotted separately. The data points in Fig. [3] 
were collected from the magnetization traces with a field 




2 4 6 8 10 12 
H(T) 



FIG. 1: (Color online) AC susceptibility, \ac, as a function 
of magnetic field H with a field sweep rate of 0.054 T/min. 
Traces are shown for T=20 mK (the highest curve), 150 mK 
and 600 mK. 

sweep-rate of 0.0068 T/min. As shown in Fig.[^a), H C 2 
saturates at 12.175 T as T approaches zero. However, 
there exists a region marked by a dashed circle at ap- 
proximately 150 mK, where a shoulder appears. The 
origin of this anomaly is unclear at this point. 

For H c i in Fig. HJb), we fit a power-law temperature 
dependence H C \(T) — H c i(0) = aT a to the data between 
1 and 260 mK, yielding an intercept of 2.053 T and an 
exponent of 1.47 ±0.1 as shown in Fig. [1] In order to 
demonstrate that the exponent is robust and indepen- 
dent of the intercept, we also performed fits to the data 
with the intercept held fixed for various intercepts be- 
tween 2.04 and 2.06 T in Fig. [3] The table in the figure 
indicates the exponent, it's error bar, and the value of % 2 
for each fit. The intercept with the lowest value of x 2 is 
-Hci(O) = 2.053 T, yielding an exponent a of 1.46 in that 
fit. Furthermore, all the fits for all the intercepts yield 
values close to 1.5. Finally we show H c i as a function 
of T 1 - 5 in Fig. 2] The inset shows the same data as a 
function of T 2 , showing that we can rule out that temper- 
ature dependence. Thus the best fit yields an exponent 
a = 1.47 ±0.1, which closely matches the expected ex- 
ponent a = 1.5 for a quantum phase transition in the 3D 
BEC universality class. Other nearby universality classes 
such as the 3D Ising (a = 2) and 2D BEC (a = 1) can 
be excluded. 

Finally, we found that DTN reached thermal equilib- 
rium at 1 mK in a few minutes which is short in com- 
parison with other solid samples that we have measured 
in the same configuration, including heavy fermion and 
two-dimensional electronic gas samples with cooling time 
constants as long as hours [26j . DTN as well as the other 
samples mentioned here were immersed in the liquid 3 He, 
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FIG. 2: (Color online) (a) The temperature dependence of the 
upper critical field H C 2 for a field sweep rate of 0.0068 T/min. 
The dashed circle shows an abnormal change in slope, (b) 
The temperature dependence of the lower critical field H c \ 
in a field sweep rate of 0.0068 T/min. The line is a fit to 
the equation H cl (T) - H cl (0) ~ T a with a = 1.47 ± 0.1 and 
Hciffi) = 2.053 T. Inset to (b): The temperature dependence 
of dXac/dH, with the critical field H c \ indicated by an arrow. 
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FIG. 3: (Color online) (a) The critical field H c \ as a function 
of temperature T. The lines are fits of the equation H c i (T) — 
Hci(0) = aT a , where H c i(0) is held fixed and a and a are 
allowed to vary. 
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FIG. 4: (Color online) The lower critical field H c i as a func- 
tion of T a , where a ~ 1.5. The dashed straight line is the 
linear fit to this trace. Inset: H c i as a function of T a , where 
a ~ 2. 



rather than being glued to a cold finger. There are several 
possible explanations for the short thermal equilibrium 
time at low temperatures. Firstly, the loose structure of 
DTN may enhance the cooling efficiency. The 3 He atoms 
can penetrate into the interlayer distance of 8.981 A. Sec- 
ondly, there exist possible exchange interactions between 
the spins of the 3 He atoms and electrons in the sample. 
Finally, a compatibility between phonon vibration modes 
in the 3 He and the sample may play a role. 

In conclusion, we have established that the field- 
induced quantum phase transition at H c \ in DTN belongs 
to the 3D BEC universality class by directly measuring 
the power-law exponent in the relation H C \(T)—H c i(§) oc 
T a down to 1 mK. This is the first example of a direct 
measurement of this exponent at temperatures far be- 
low the energy scales for antiferromagnetic coupling in a 
magnetic insulator. 

This work was carried out at the High B/T facility of 
the National High Magnetic Field Laboratory supported 
by the National Science Foundation Cooperative Agree- 
ment No. DMR 0654118 and the State of Florida. 
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